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Asteroseismic determination of helium abundance in 
stellar envelopes 
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ABSTRACT 

Intermediate degree modes of the solar oscillations have previously been used to de- 
termine the solar helium abundance to a high degree of precision. However, we cannot 
expect to observe such modes in other stars. In this work we investigate whether low 
degree modes that should be available from space-based asteroseismology missions 
can be used to determine the helium abundance, Y, in stellar envelopes with sufficient 
precision. We find that the oscillatory signal in the frequencies caused by the depres- 
sion in Ti in the second helium ionisation zone can be used to determine the envelope 
helium abundance of low mass main sequence stars. For frequency errors of 1 part in 
10 4 , we expect errors ay in the estimated helium abundance to range from 0.03 for 
0.8M Q stars to 0.01 for 1.2M© stars. The task is more complicated in evolved stars, 
such as subgiants, but is still feasible if the relative errors in the frequencies are less 
than 10~ 4 . 
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1 INTRODUCTION 

The helium abundance of stars is often used to extrap- 
olate back to obtain estimates of the primordial helium 
abundance which is an important test of cosmological mod- 
els. The helium content of the oldest stars cannot how- 
ever, be measured directly — these are low mass stars, 
and their photospheres are not hot enough to excite he- 
lium lines that can be used to determine the helium abun- 
dance spectroscopically. Thus the helium content of low- 
mass stars has to be derived from the evolution of heav- 
ier elements via stellar models. The evolution of stars how- 
ever depend crucially on the helium content itself, mak- 
ing the knowledge of the helium abundance in stars very 
important. The estimated ages of globular clusters, which 
provides an important constraint on cosmological models, 



diffusion in the envelope (Delivannis, Dcmaraue & Kawalcr 


1990; 


Proffitt & VandenBcrg 


1991; Delivannis & Dcmaraue 


1991; 


Chabovcr et al.lll992|). 



The helium abundance of the solar envelope h as been 
successfully determined using helioseismic d ata iGoughl 
Il984l; iDappen et a.ljTl99lt iBasu fc Antialll995l etc.) These 
seismic determinations were based on the fact that the ion- 
isation of hydrogen and helium causes a distinct, localised 
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depression in the adiabatic index, Y\, in the near-surface lay- 
ers of the Sun. The first helium ionisation zone is not very 
useful in making seismic inferences — it overlaps with the 
hydrogen ionisation zone and is located very close to the 
stellar surface where there are significant uncertainties in 
seismic inversions. The second helium ionisation zone which 
is located below the highly superadiabatic layer of the con- 
vection zone is not particularly sensitive to surface effects 
and hence is useful in determining the helium abundance. 
The depression of Pi in this region can be directly related 
to the abundance of helium; this depression increases with 
increase in the helium abundance (see Fig. 0. The depres- 
sion in Pi causes a localised depression in the derivative of 
the sound speed in the region, and hence affects all acoustic 
modes that travel through that region. Thus using helioseis- 
mic inversion techniques it is possible to determine the sound 
speed in the solar interior which in turn could be calibrated 
to obtain the helium abundance in the solar envelope. All 
studies to determine the helium abundance in the Sun used 
intermediate degree modes to measure the sound speed in 
the Hell ionisation zone. 
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Figure 1. The adiabatic index, Ti, as a function of fractional 
radius for 1.2 Mq models. The upper two panels are models with 
different Y evolved to the age indicated in the panels. The bottom 
panel shows Ti for 1.2 Mq models evolved till the surface is at 
a radius of 1.15-Rq. The lines of different types represent models 
with different helium abundances as marked in each panel. 

sufficient accuracy to enable a more detailed study of stellar 
structure than has been possible so far. Ground based 
observations have already started measuring the frequencies 
of many stars desp i te the difficulties of such measurements 
feouchv fc Carrier! l2002t) . Unfortunately, none of these 
observations can or will be able to observe the intermediate 
degree modes of oscillation that have been so useful in solar 
studies. These missions can at most, hope to determine 
the frequencies of modes with degree £ ^ 3. Therefore, it 
is important to devise a practical way for determining the 
helium content of stellar envelopes using only low degree 
modes. 

Any localised feature in the sound speed inside a star, 
such as that caused by the change in the temperature gra- 
dient at the base of the convection zone, introduces an os- 
cillatory term in frequencies as a function of radial order n, 
that is proportional to 



sin(2r TO w n ,« + 4>), 




(e.g., iGouehl Il990l) where r m is the acoustic depth of 
the localised feature, c the speed of sound, r m the ra- 
dial distance where the feature is located; and u> n ,e the 



angular frequency of a mode with radial order n and 
degree I, <f> is a phase. This oscillatory signature has 
been extensively studied for the Sun in order to deter- 
mine the extent of overshoot below the solar convection 
zone l|M onteiro 1 Christensen-D alseaard ~ Thompson! Il994l : 
iBasu. A^nda^^Narasmihaf^^9^ . The Hell ionisation zone 
too gives rise to a similar oscillatory signal, which can be 
extract ed from frequencies of low deg ree modes (see for ex- 
ample, iRoxburgh fc Vorontsovl 1200111 . The present work is 
motivated by the possibility of using this signal in low de- 
gree seismic data from distant stars to study the properties 
of the Hell ionisation zone. 

The depression in T\ in the Hell ionisation zone is not 
a very sharp feature and has a finite width. However, one 
can show that it still gives rise to an oscillatory signal in the 
frequencies and that the amplitude of the oscillatory term 
caused by the depression in Yj depends on the amount of he- 
lium present at the location of ionisation. In principle, if this 
amplitude is measured, it can be compared against models 
to estimate the helium abund ance in the stellar envelope. 
iMonteiro fc Thompson! jl99ct) have studied the oscillatory 
signal arising from the Hell ionisation zone in stellar mod- 
els. They have shown that because of the finite width of 
this zone, the amplitude of oscillatory signal in frequencies 
is modulated by an additional factor 

sin 2 M)/M)> (2) 

where j3 is half the acoustic thickness of the dip in I?i. Since 
(3 is rather small, of the order of 100 s, this factor pro- 
vides a slow modul ation in amplitude of the oscillatory term. 
iMiglio et al] i200St) have also examined the oscillatory signal 
due to Hell ionisation zone and pointed out the possibility 
of u sing it to determine the helium abundance. 

IBasu. Antia fc Demarauel i2003h examined several stel- 
lar models to study the possibility of determining helium 
abundance from the amplitude of the oscillatory signal. 
They calculated the amplitudes using the second differences 
of the frequencies. In this work we extend that study to a 
much larger sample of models to explore possible systematic 
errors in the measurement of helium abundance using only 
low degree modes. We also look at models with substantially 
different abundances, such as stars with Population II abun- 
dances. We also look at stars beyond the main sequence. The 
rest of the paper is organised as follows: we describe the 
measurement method in Section and give details of the 
models we have used in this study in Section [3] the results 
(including those for solar low degree data) are presented in 
Section 01 and we discuss our conclusions in Section |S] 



2 THE TECHNIQUE 

The amplitude of the oscillatory signal due to the depression 
in Ti in the Hell ionisation zone is very small and it may 
be desirable to amplify it for a proper measurement. We 
amplify the signal by taking the second difference of the 
frequencies, 

S 2 V n ,l = Vn+\,t — 1l>n,l + Vn-l,l (3) 

for modes with the same degree. We use modes of degrees 
£ = 0-3. In addition to amplifying the oscillatory signal, the 
process of taking second differences suppresses the dominant 
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smooth trend of the frequencies to a large extent, making 
it easier to fit and measure the oscillatory term. Th e sec- 
ond d ifferences are then fitted to the form used by iBasul 
Jl997f) . but without the degree-dependent terms (which are 
not relevant for low degrees), i.e.; 

8 2 u = fai + a 2 u + 

+ (bl + ^|^) SUl(47r^THo + 0Hc) 

+ ( Cl + i^) sin (^ 7TUTCZ + 0cz), (4) 

where the first term with coefficients 01,02,(13 defines the 
smooth part of the second difference, the second term is the 
oscillatory signal from the Hell ionisation zone, and the last 
term is the oscillatory signal from the base of the convection 
zone. In Eq.0] the terms 7~He and tcz are the acoustic depths 
of the Hell ionisation zone and the base of the convection 
zone, respectively. The terms <f>Ke and <j>CZ are the phases of 
the corresponding oscillatory signals. The parameters 01, 02, 
03, bi, 6 2 , t Hc , 0Ho, ci, c 2 , rcz and (pcz are determined by 
least-squares fits to the second differences of the frequencies. 

As mentioned earlier, the Hell ionisation zone actually 
has a small but finite width and the first derivative of the 
squared sound speed (after some scaling) can be better ap- 
proximated by a triangular function rather than a disconti- 
nuity. To account for the extra modulation in the oscillatory 
term due to the finite width of the ionisation zone (see Eq.EJ 
we have also fitt ed the oscillation with the mo dulating factor 
sm 2 (LuP)/(ujP) iMonteiro fc Thompsonlll99^ in the ampli- 
tude keeping /3 as a free parameter: 

5 2 v — (ai + aiv + azv 2 } 

( & lS in 2 (27n//3)\ . fA 
+ I I Sin(47ri/THe + </>He) 

+ ( Cl + ~ + ~l) sin (47r^-rcz + 4>cz), (5) 

For most of the cases considered in this paper, f3 is found 
to be of the order of 100 s and hence, for the range of fre- 
quencies included in the fits, the term sm 2 {2-Rv(3) can be 
approximated reasonably well by a quadratic term as as- 
sumed in Eq. (J1J. Thus although the fitting form in Eq. (|1J 
had been empirically determined to fit the signature from 
the base of the solar convection zone using intermediate de- 
gree solar modes, it is expected to fit the signature from 
Hell ionisation zone in low degree modes quite well. 

To study the sensitivity of fitted amplitude and acoustic 
depth to the fitting form we have tried one more variation 
where the oscillatory term due to Hell ionisation zone in 
Eq. JSJ is modified to 

/ 62 &3 \ 

( bi + — + J sin(47r^r Hc + 0h c ) • (6) 

This form is obtained by approximating the factor 
sin (2irv0)/{yj3) in Eq. @ by a quadratic in \/v. Since j3 is 
small such an approximation is reasonably good. This form 
differs from Eq. Q through the addition of a linear term in 
1/v in amplitudes. 

The fits of the three forms to frequencies generated by a 
I.2M0, 1.2i?0 model are shown in Fig. El Each fit is a com- 
posite function representing the smooth trend in the second 
differences, a dominant slow oscillation due to the second 
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Figure 2. A sample of the fits to the second differences of the 
scaled frequencies. The points are the 'data', and the lines the 
fits to the points. The examples shown are for a 1.2Mq model 
evolved to a radius of 1 . 2Rq . 

helium ionisation zone, and a higher frequency oscillation 
arising from the base of the convection zone. Each panel in 
the figure shows two fits, one for each of the two extreme val- 
ues of the initial helium abundance Yo that we have chosen 
for this study. As can be seen in the figure, the amplitude 
of the slower oscillation clearly increases with the amount 
of helium present in the envelope. We aim to calibrate this 
increase in the amplitude of the oscillatory signal from the 
Hell ionisation zone against the helium content. We do such 
fits for each of the stellar models. Since the amplitude of the 
signal is a function of frequency, we use the fits to obtain 
the mean amplitude over the fitting interval. There is not 
much difference in the fits obtained using the three forms 
given by Eqs. JlJJSJl, as borne out by Fig. This ensures 
that the amplitudes are not too sensitive to the particular 
fitting form chosen. 

The frequencies of stellar oscillations roughly scale as 
y/p, where p tx M/R 3 is the mean density of the star. Since 
the models have different mean densities, it is convenient to 
scale the frequencies by the factor 

/ = Vp/P^ (7) 

before taking the differences and fitting the signal. How- 
ever, this requires a knowledge of M/R 3 for the stars to be 
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studied, which may not be known for all stars that may be 
targets of asteroseismic study. This factor can actually be 
estimated from the seismic data using the large frequency 
separation, 

= v n +i,i — fn,i (8) 

(see e.g., iMazumdar fc Antial l200ll : IChristensen-Dalsgaardl 
2003). We find that the averaged large frequency separation 
scaled by the factor in Eq. £J i s around 135 /iHz for all the 
main-sequence stellar models in the range of masses consid- 
ered in this study. We therefore, scale the frequencies using 
the average frequency separation to avoid the use of the po- 
tentially unknown factor M/R 3 . We have chosen a frequency 
range of 1.5 mHz to 3.0 mHz after scaling by A/A©, where 
A is the mean large frequency separation and Aq = 135 ^tHz 
is the value for the Sun. The choice of the frequency range 
is dictated by practical reasons. At higher frequencies the 
oscillatory signal due to the Hell ionisation zone becomes 
small. Furthermore, the observational error is expected to 
increase in higher frequency modes, as is the case for solar 
data. At lower frequencies there is significant departure from 
the fitting forms assumed in this work since the low n modes 
can not be approximated well by the asymptotic for mulae 
used in obtaining the oscillatory form l|Monteirolll99ft) . 

It can be shown that the process of taking the second 
differences magnifies the amplitude of the oscillatory func- 
tions by a factor 

4sin 2 (27rrA) (9) 

where A is the mean large frequency separation. In order to 
obtain the real amplitude of the oscillatory signal, i.e., the 
amplitude of the oscillations in the frequencies, we divide the 
amplitude obtained for the second differences by this fac- 
tor. This conversion back to frequencies is found to remove 
some variation in the am plitude between different models 
iMazumdar fc Antiall200lT) . All the amplitudes presented in 
this work are the amplitudes of the oscillatory function in 
the frequencies that we obtain by converting the fitted am- 
plitudes in the second differences using the conversion factor 
in Eq. @. We find that for the Hell ionisation zone, this fac- 
tor is not very different from unity and hence the amplitudes 
are similar for both frequencies and the second differences 
of frequencies. This magnification is however quite large for 
the signal from the base of the convection zone, making it 
easier to fit the two oscillatory terms simultaneously in the 
second differences. 

By fitting the oscillatory term in the second differences, 
it is possible to determine the characteristics of the Hell ion- 
isation zone, like its acoustic depth and the thickness and 
amplitude of the depression in Ti from the fitted parameters. 
One would expect the amplitude of the oscillatory signa- 
ture to increase with helium abundance and he nce it can be 
calibr ated using stellar models with known Y. iMiglio et alJ 
( 2003) have suggested that the area of the bump in Ti should 
be a measure of helium abundance. The product of ampli- 
tude with half-width (3 provides a measure of this area. We 
find that it is difficult to estimate j3 reliably by fitting the 
oscillatory term. We find a large scatter in the results if the 
product of the amplitude with j3 is used, which makes it 
difficult to determine Y and hence we do not use it in our 
study. Instead, we try to find a relation directly between the 
amplitude of oscillatory term and the helium abundance. 



3 STELLAR MODELS 

We use the technique described in the previous section to 
study the properties of the oscillatory signal in the frequen- 
cies of low degree modes and its relationship with helium 
abundance in the stellar envelope using a large number of 
stellar models. We restrict ourselves to stars with masses 
close to the mass of the Sun, namely, the range 0.8-1.4 M®. 
The upper limit of the mass range is set by the fact that it 
is difficult to extract the oscillatory signal in more massive 
stars. This happens because the acoustic depths of the base 
of the convection zone and that of the Hell ionisation zone 
are very similar in higher mass stars, and the two oscillatory 
signals, one from the convection zone base and the other 
from the helium ionisation zone, interfere with each other. 
The lower mass limit is dictated by the fact that the am- 
plitude of the oscillatory signal reduces with mass, though 
there is no particular difficulty in determining the signal as 
we are using exact model frequencies. In actual practice the 
lower limit on the stellar mass that can be studied will be 
determined by the precision to which the frequencies can 
be determined from the future asteroseismic missions. Since 
observed frequencies will have uncertainties associated with 
them, we have examined the effect of frequency errors on 
the fits in Section 2] We study stellar models evolved to dif- 
ferent ages on the main sequence for masses 0.8Mq, IMq, 
1.2Mq and 1.4M©. We have also constructed a few stellar 
models which have evolved beyond the main sequence and 
are on the subgiant branch. 

The stellar models were constructed using YREC, the 
Yale R otating Evolution Code in its non-rotating configu- 
ration iGuenth er e^dj^)92|). These models use the OPAL 
equa tion of state jRogers fc Navfonovll2002T) . OPAL opaci- 
ties ijlglesias fc Ro ger! J^^)7"k )w temperature opacities of 
lAlexand er fc Ferguson ljl994h and nucle ar reaction rates as 
used bv iBahcall fc Pinsonneaultl il992T) . The models take 
into account diffus ion of helium and heavy elem ents, using 
the prescription of lThoul. Bahcall fc Loebl dl994h . 

Most of the models have been constructed with Pop- 
ulation I abundances. These models have an initial heavy 
element abundance Zo = 0.022 and were constructed with a 
mixing length parameter a — 2.1. The initial heavy element 
abundance and the mixing length are the same as those used 
to construct ca librated standard solar models using YREC 
iWinnick et all 2002). The value of the initial helium abun- 
dance Yo is varied: we use 5 values of Yq for the Population I 
models, Y = 0.24, 0.26, 0.28, 0.30 and 0.32. Diffusion and 
gravitational settling of helium causes the helium abundance 
in the stellar envelope to change with time. Therefore, the 
helium content in the envelope of a star still evolving on 
the main sequence will decrease progressively with age from 
its initial (or zero-age main sequence) helium abundance Yo. 
Beyond the main sequence the dredge-up of helium begins 
to increase the helium abundance in the envelope. Thus the 
amplitude of the Hell signal will not reflect the initial abun- 
dance Yo of the models; it will depend on the abundance 
of helium in the stellar envelope at the stage we examine 
it. As a result we always use the actual, rather than the 
initial, helium abundance in the envelope to compare and 
calibrate the amplitude of oscillatory signal. To investigate 
systematic errors caused by model parameters, we have also 
constructed models with a slightly different value of the ini- 
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Figure 3. The scaled acoustic depth th c and the mean amplitude of the oscillatory part in frequency due to the Hell ionisation zone in 
the scaled frequency range of 1.5-3.0 mHz are shown as a function of the helium abundance in the envelope. The open symbols represent 
0.8AfQ models, skeletal symbols show XMq, starred for 1.2Mq and filled symbols show 1.4Mq models. For each model the results for 
the three independent fits using Eqs. J3IE1EJ are shown by 3, 4 and 5 sided symbols respectively. 



tial heavy element abundance (Zo = 0.018) and some models 
with a slightly different value of the mixing length parameter 
(a = 1.7). 

To test whether the method outlined in this paper also 
applies to stars of very different heavy element abundances, 
we have two other sets of models. One set with initial abun- 
dance Zo of 0.007 (the abundance of stars like r Ceti) and 
another set with Population II abundance of Zo — 0.001. 
Models for these two sets were constructed with initial Y of 
0.22, 0.24, 0.26, 0.28 and 0.30. For each of the stellar mod- 
els we calculate the frequencies of low degree modes and 
attempt to extract the oscillatory term as explained in the 
previous section. 

From our experience of working with solar data, we ex- 
pect that one of the important sources of uncertainty is the 
equation of state (EOS) which determines the width and the 
depth of the dip in Fi . The role of the equati on of state has 
been extensively studied f or the solar case iAntia fc Basul 
11994 iRichard et aljll998l lBasulll99St) by looking at mod- 
els constructed with OPAL, t he Mihalas, Hummer fc 
D&ppen (henceforth MHD1 EOS |Hummer fc MihalaJl988t 
iMihalas. Dappen fc Hummerlll98St iDappen et alJll988T) as 
well as other simpler EOS. lMonteho*fcTnomra5QiTll998h 
have studied the influence of equation of state in stellar 
case too and find that the half-width j3 of the Hell ioni- 
sation zone is determined by the EOS. Helioseismic studies 
have shown that in the helium ionisation zone, the OPAL 
EOS is a better re p resentation of stellar EOS (see e.g., 
Basu fc Antial Il995t iBasu fc Christensen-Dalsgaardl Il997t 
BastT^DaOTemfc^a^bnovf^^99l) and hence, in this work 
we restrict ourselves to models with OPAL EOS, though we 
do study the effect of the EOS using a few IM@ models con- 
structed with the MHD equation of state. This should give 
a reasonable estimate of the uncertainty associated with the 
EOS. 



4 RESULTS 

4.1 Population I Main Sequence Models 

Fig- El shows the acoustic depth and the average amplitude 
of the oscillatory signal due to the Hell ionisation zone for 
the standard set of main sequence models used in the study. 
It can be seen that the acoustic depth th c is within a factor 
of two for all models, while the amplitude varies over a large 
range. It should be noted that all values of the acoustic depth 
and the amplitude listed in this paper are those obtained af- 
ter scaling by the frequency separation. The actual values of 
7~He may show significant variation due to variation in M/R 3 
or equivalently in A. On average, the amplitude appears to 
increase with Y, the helium abundance in the stellar enve- 
lope, but for a given helium abundance, the amplitude also 
increases with stellar mass. This figure combines the results 
obtained using different fitting forms and it can be seen that 
there is a good agreement between different results, though 
for 1.4Mq models the difference tends to become large as 
the fits are not very good for many of these models. The 
mean difference in amplitudes obtained using two different 
fitting forms is found to be about 0.02 /iHz for stellar models 
with masses O.8M0, while for IMq and 1.2M© models it is 
about 0.06 /iHz. As we shall see later, these differences are 
smaller than the errors that may be expected to arise from 
uncertainties in observed frequencies. For I.4M0 models the 
mean difference in amplitudes between two fitting forms is 
about 0.18 /iHz, which is comparable to the errors expected 
to be caused by uncertainties in the observed frequencies. To 
show the variation of amplitude with Y more clearly, Fig. 2] 
shows the amplitude as a function of Y for each mass sepa- 
rately. There is considerable scatter around the mean trend 
due to factors like age, and hence the evolutionary state. 

The dependence of the amplitude on factors other than 
the helium abundance complicates the determination of the 



6 S. Basu et al. 




Figure 4. The mean amplitude of the oscillatory part in frequency due to the Hell ionisation zone is shown as a function of the helium 
abundance in the envelope. Each panel shows the result for models with a fixed mass. The triangles, squares and pentagons show the 
results using fitting functions JFJ and J1JJ respectively. 
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Figure 5. The mean amplitude of the oscillatory part in frequency due to the Hell ionisation zone is shown as a function of the helium 
abundance in the envelope for models with fixed mass and radius. The different masses and radii are marked in each panel. The solid, 
dotted and dashed lines show the results for fits using Eqs. HI . 151 and JHJ respectively. For sake of clarity, estimated errors are shown 
on only one set of points. 



helium abundance of a given star. Our task is made easier 
if conventional observations combined with modelling give 
us a reasonable idea of the radius, R and mass M of the 
star. If only the radius is known, we can estimate M/R 3 
from the frequency separation and combining that with the 
known radius will give an estimate of mass. Similarly, if only 
the stellar mass is known independently, we can use the fre- 
quency separation to estimate the stellar radius. This would 
enable us to make calibration models with a given radius 
and mass. A plot of the mean amplitude as a function of Y 
for models with given masses and radii is shown in Fig. |S| 
In this case all we would need to do to determine Y is to 
determine the mean amplitude from the observations and 
interpolate. Models of the same mass but different initial 
helium abundances will have different ages when they have 
the same radius, as we shall see in section [Ql this does not 
affect the calibration process. 

Fig. also shows an estimate of the uncertainties in Y 
expected as a result of observational errors. Shown in the fig- 
ure are the error bars on a few representative points, which 
are obtained assuming a constant relative error of 1 part in 
10 4 in the frequencies. This level of error is not unrealis- 
tic given that even ground-based observations have almos t 
reached this level of precision (see iBouchv fc Carriej Eo02) . 
These error bars have been calculated using a Monte-Carlo 
simulation where random errors with specified standard de- 
viation were added to the exact frequencies before fitting 



the data. The results of the fits to the error-added frequen- 
cies were then used to estimate the variance in each fitted 
parameter to calculate the expected error. The precision to 
which we can use the amplitude to determine Y increases 
with increase in mass and age. From the error bars it ap- 
pears that the expected error in Y would range from less 
than 0.01 for 1.2M models to about 0.03 at 0.8M Q . To get 
better precision for 0.8Mg stars we will need to measure the 
frequencies to a relative accuracy of better than 10 -4 . The 
error in the scaled acoustic depth r of the Hell ionisation 
zone is about 20-40 s except for 0.8Mq models where it is 
about 100 s. The above results were obtained assuming we 
had modes with degrees of 0,1,2 and 3. The uncertainty in 
the Y determination increases if I = 3 modes are unavail- 
able, though the increase is very marginal and will not affect 
any of the results. Obviously, there would be additional er- 
rors due to uncertainties in our knowledge of stellar mass 
and radius. This error can be estimated by considering all 
models within the permitted range of mass and radius. 

It should be noted that the points shown in figures have 
been obtained using the exact model frequencies, only the 
extent of the error bars were obtained from the Monte-Carlo 
simulations. From simulations it is found that when random 
errors are added to the frequencies, the estimated ampli- 
tudes increase by an amount comparable to the estimated 
errors. T his effect ha s been seen earlier for helioseismic data 
(see e.g.. lBasulll997i) . Therefore for the purposes of calibra- 
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Figure 6. The mean amplitude of the oscillatory part in fre- 
quency due to the Hell ionisation zone is shown as a function of 
the helium abundance in the envelope for lMg stars that have 
a radius of IRq. The results obtained for models with a slightly 
different initial heavy element abundance Zo and mixing length 
parameter a are compared. Representative errors in the ampli- 
tudes are indicated on one set of points. All results were obtained 
by using the function in Eq. JSJ to fit the data. 
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Figure 7. The mean amplitude of the oscillatory part in fre- 
quency due to the Hell ionisation zone is shown as a function 
of the helium abundance in the envelope. Each panel shows the 
results for a given bin in the ratio of the large frequency sepa- 
ration A/Aq. For clarity only results obtained using Eq. JSJ arc- 
shown. Open squares, crosses, stars and filled squares displays 
stellar models with masses O.8M , IMq, 1.2Mq and 1.4M Q re- 
spectively. 



tion once real data are available, the calibration curve should 
be obtained after adding the random errors to the frequen- 
cies of the calibration models. We have not done that in this 
work as in that case the figures will depend somewhat on 
the estimated errors in frequencies. 

To study the influence of the initial heavy element abun- 
dance, Zo, and of the mixing length parameter, a, on the 
amplitude of oscillatory signal we have constructed stellar 
models with different values Zo and a for stars with the 
same masses and radii shown in Fig. Figure HJ shows the 
results for IMq stars that have the same radius (17?©) but 
have different Zo and a. It can be seen that the variation in 
amplitude due to reasonable uncertainties in Zo and a are 
smaller than the variation due to Y , and also smaller than 
the error estimates for a relative error of 10 -4 in frequencies. 

In the above discussion we have assumed that either 
the mass or the radius of the star is known independently. 
This is a reasonable assumption since most of the targets 
of asteroseismic missions are likely to be nearby stars for 
which at least, one of these quantities is known. However, if 
for some star neither the mass nor the radius is known, then 
from the frequency separation we can only estimate the ratio 
M/R 3 . To examine the implication in this case, we show in 
Fig. [7] the amplitude as a function of Y for all models in a 
given bin of the mean large frequency separation A. It can be 
seen that there is some scatter due to variation in mass and 
radius, but in most cases the trend of increasing amplitude 
with Y is very clear and can be calibrated to calculate the 
helium abundance. Of course, the systematic errors will be 
larger in this case. 



4.2 Results for the Sun 

In order to test this technique we can try to estimate the 
helium abundance in the solar envelope using only the 
low degree modes. We use the frequencies obtained from 
the first 360 days of observations by the MDI instr ument 
on board the SOHO spacecraft llSchou et al]ll99Sft . The 
observed frequencies with degrees of £ 3 and frequencies 
in the range 1.5 ^ v 3.0 mHz were fitted to obtain the 
mean amplitude of the oscillatory signal due to the Hell 
ionisation zone, and the result is shown by the dotted 
horizontal line in Fig. |H| To infer the helium abundance in 
the Sun we have used IMq, IRq models for the calibration. 
It should be noted that while these models have the same 
mass and radius as the Sun, these calibration models are 
not solar models, the models have different ages and lumi- 
nosities. To the frequencies of each model we added fifty 
realisations of the observed errors, and for each realisation 
we determined the amplitude of the oscillations. We have 
plotted the mean amplitude for each model as a function 
of the envelope helium abundance of the model in Fig. |H| 
All these fits use the same set of modes that were used 
to extract the oscillatory signal from the observed data. 
To study the sensitivity of results to uncertainties in the 
EOS, the figure also shows the results of a set of models 
constructed with the MHD equation of state. The result of 
the observations is shown as the horizontal line and we can 
see that we obtain a helium abundance of 0.239 ±0.005 with 
the MHD models and 0.246 ± 0.006 with the OPAL models. 
Similar results have been obtained using frequencies from 
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Figure 8. The mean amplitude of the oscillatory part in fre- 
quency due to the Hell ionisation zone is shown as a function of 
the helium abundance in the envelope for models with M = IMq 
and R = IRq . The horizontal line shows the results obtained us- 
ing observed frequencies from MDI instrument for the Sun. The 
two lines show the results using OPAL and MHD EOS as marked 
in the figure. These results are obtained using the function in 
Eq. to fit the data. 



GOLF. These results agree with the results obtained using 
interm e diate degree solar p-modes (e.g.. lAntia fc Basu l 
[r99^^_|Pe£ez Ifemandez fcT^risto 1 1994 

basu fc AntisJ 1199.4 lBasulll99Sl: iRichard et al]ll99Sfr . This 
exercise therefore demonstrates that it is indeed possible to 
determine helium abundance using only low degree modes. 



4.3 Stellar models beyond the main sequence 

All the stellar models that we have considered so far are on 
the main sequence. We have constructed a few models which 
have evolved beyond the main sequence to see how we might 
determine the envelope helium abundance for evolved stars. 
For these stars, many nonradial modes show a mixed charac- 
ter where the same mode behaves like a gravity mode in the 
core and like an acoustic mode in the envelope. The analy- 
sis that leads to the calculation of the oscillatory signal in 
frequencies is not valid for such modes. In fact, even cal- 
culating the large frequency separation from the calculated 
frequencies is somewhat difficult for such stars. However, 
radial modes are assured to be purely acoustic in nature 
throughout the star, and hence, for such stars we restrict 
our analysis to radial modes [i — 0) only. Thus the number 
of modes available for fitting is small and we do not expect 
to be able to calculate all parameters in Eqs. JIJ, 10 or © 
reliably. Since the amplitude of the oscillatory signal due to 
the base of the convection zone is generally much smaller 
than that due to the Hell ionisation zone, we only fit the 
signal due to the Hell ionisation zone. Sample fits are shown 
in Fig. The fits are reasonably good, though not as good 
as those in Fig. [2] hence we expect somewhat larger system- 
atic errors. Since the amplitude of the Hell signal increases 
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Figure 9. A sample of the fits to the second differences of the 
scaled frequencies using only 1 = modes for models evolved 
beyond the main sequence. The points are the 'data', and the 
lines the fits to the points. The examples shown are for 1.2Mq 
model evolved to a radius of 1.8-Rq. 



with reduction in frequency, we have extended the fitting 
interval by including modes with frequencies 1-3 mHz. The 
enhanced range increases the number of modes available as 
well as the mean amplitude of the signal, thus providing bet- 
ter characterisation of the signal. When compared to main 
sequence stars, we expect larger uncertainties arising from 
errors in frequency because of the fewer number of modes 
in these stars. Nevertheless, Monte-Carlo simulations with 
comparable error in frequencies show that the estimated er- 
ror in mean amplitude for evolved stars is not too much 
larger than that for main-sequence stars. 

Figure [TU1 shows the amplitude of the signal as a func- 
tion of Y for models with a fixed mass and radius. Error 
bars corresponding to frequency errors of 1 part in 10 4 are 
also shown for a few typical models. As expected the am- 
plitude increases with Y and it is possible to determine Y 
from known amplitude. We find again that these amplitudes 
are not particularly sensitive to small variations in Zo and 
a. However, in this case the differences due to Zo and a are 
larger than those for main sequence models. As can be seen, 
despite the errors, we should be able to determine the he- 
lium abundance in the envelope of these stars to a precision 
of about 0.01-0.02. 
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Figure 10. The mean amplitude of the oscillatory part in frequency due to the Hell ionisation zone is shown as a function of the helium 
abundance in the envelope for evolved stellar models with a fixed mass and radius. The solid, dotted and dashed lines show the results 
for fits using Eqs. J5J and JSJ respectively. The left panel shows the results for a few different masses and radii as marked in the 
figure, while the right panel shows the results for M = IMq, R = 1.75Bg models with different Zo and a. For sake of clarity, the error 
estimates are shown on only one set of points. 



4.4 Results for models with low heavy element 
abundances 



Stars with very low heavy element abundances evolve 
quite differently from stars with near-solar heavy element 
abundances. This is because the convection zone of low 
Z stars is shallower a t the same effective tem perature 
JPelivannis et al.l Il99(j iProffitt fc Michaudl fl99lh . We do 
not expect this to affect the helium signature; nevertheless, 
we confirm this using models of low Zo . We have constructed 
two sets of models, models with Zo = 0.007 and models with 
Zo — 0.001. We have both main sequence and subgiant mod- 
els with Zo = 0.007, but only subgiants with Zo = 0.001 
since Population II models evolve very quickly compared 
with stars of the same mass with higher Z and hence we 
do not expect to see many main-sequence stars. As with the 
Zo = 0.022 models, we use £ = 0,1,2 and 3 modes for the 
main sequence stars and only £ = modes for the subgiants. 

Fig. 1111 shows the amplitude as a function of Y for 
Zo = 0.007 models. Again we can see that the amplitude 
increases with Y . The amplitudes are however slightly dif- 
ferent from those of Zo = 0.022 models, and hence for such 
large differences in Zo it will be better to make different cal- 
ibration models. Fig. 1111 also shows the results for Popula- 
tion II models. Again the amplitude increases with Y. The 
amplitude for these evolved models are much larger than 
those for main sequence models studied in Sect. 14.11 Nev- 
ertheless, to get good results for Y, we would need lower 
errors in frequencies. We expect to be able to observe only 
lower mass Population II stars, and this implies that the 
expected amplitude of the oscillation signal will be small, 
making errors larger. 



5 CONCLUSIONS 

We have presented results of our investigation into how low 
degree oscillation frequencies may be used to determine the 
helium abundance in stars with solar-type oscillations. We 
have demonstrated that in principle, it is possible to use 
this technique to estimate the helium abundance in stellar 
envelopes using the frequencies of low degree modes of os- 
cillation. 

We find that the oscillatory signal in the frequencies 
due to the helium ionisation zone can be used to determine 
the helium abundance of a low-mass main-sequence star, 
provided either the radius or the mass is known indepen- 
dently. The precision to which we may be able to determine 
the helium abundance increases with increase in mass. Us- 
ing reasonable error estimates it appears that a precision 
of 0.01-0.02 in Y is possible in most cases. The amplitude 
of the oscillatory term is not particularly sensitive to small 
variations of the initial heavy element abundance, Zo and 
mixing length parameter, a. For stars with masses of about 
1.4Mq and larger it is difficult to fit the oscillatory signal 
reliably and hence the uncertainties are larger. This diffi- 
culty is most probably due to the fact that these stars have 
shallow convection zones and there is not much difference in 
the acoustic depths of the Hell ionisation zone and that of 
the base of the convection zone, thus leading to interference 
between the two terms. For low mass stars the amplitude 
of the oscillatory term is small leading to large errors in de- 
termining the helium abundance. Thus the best results are 
obtained for stars with masses close to and larger than the 
solar mass. Using solar data, we have shown that we can 
indeed determine the envelope helium abundance using low 
degree modes. 

Even if the mass and radius for a star is not known 
independently, we can use the large frequency separation 




Figure 11. The mean amplitude of the oscillatory part in frequency due to the Hell ionisation zone is shown as a function of the helium 
abundance in the envelope for stellar models with fixed mass and radius, constructed with low heavy element abundance. This figure 
shows the results for stellar models with initial low heavy element abundance, Zq = 0.007 (left panel) and Z§ = 0.001 (right panel). The 
solid, dotted and dashed lines show the results for fits using Eqs. and JHJ respectively. For sake of clarity, the error estimates are 

shown on only one set of points. 



A to estimate M/R 3 and in that case also a reasonably 
tight correlation is found between Y and amplitude of the 
oscillatory signal, thus allowing us to determine Y, though 
with larger errors. 

For stellar models that are evolved off the main se- 
quence it is not possible to use nonradial modes, however by 
using only radial modes it is possible to fit the oscillatory 
term arising from the Hell ionisation zone, which can then 
be used to measure the helium abundance. A more detailed 
study of stars in different post main sequence phases is re- 
quired to understand the oscillatory signal from such stars. 
Models with much lower heavy element abundances behave 
in a manner similar to those with near-solar heavy element 
abundances as far as the helium signature is concerned. 

Although most of the present work uses exact frequen- 
cies from stellar models, we have successfully applied this 
method to solar low degree data. We have also estimated the 
expected errors assuming frequency errors of one part in 10 4 , 
a level that is expected to be achieved by future space- and 
ground-based observations. This opens up the possibility of 
using asteroseismic data to obtain the helium abundance in 
stellar envelopes, leading to a better understanding of the 
evolution of stellar populations in our galaxy. 
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